Spin-triplet superconductivity due to antiferromagnetic spin-fluctuation in Sr 2 Ru04 
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A mechanism leading to the spin-triplet superconductivity 
is proposed based on the antiferromagnetic spin fluctuation. 
The effects of anisotropy in spin fluctuation on the Cooper 
pairing and on the direction of d vector are examined in the 
one-band Hubbard model with RPA approximation. The gap 
equations for the anisotropic case are derived and applied to 
Sr2Ru04. It is found that a nesting property of the Fermi 
surface together with the anisotropy leads to the triplet su- 
perconductivity with the d = z(sink x ± isink y ), which is 
consistent with experiments. 
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Since the discovery of superconducting phase in 
Si'2Ru04 much effort has been paid for understand- 
ing its exotic properties. Among several interesting na- 
tures, the most fascinating one is that it is a spin- 
triplet superconductor confirmed by NMR experiment 
While most superconductors found during several 
decades are singlet, the only exceptions were 3 He and 
UPt3. Therefore the fact that the triplet pairing is re- 
alized in Sr2RuC>4 has attracted much attention. While 
UPt3, the second example of spin-triplet superconduc- 
tor, has a complicated electronic structure, Sr2Ru04 has 
a rather simple electronic state Qj. Thus clarifying the 
microscopic mechanism of superconductivity in Sr2Ru04 
is very important for understanding the triplet supercon- 
ductors in general. 

In 3 He, Cooper pairs are formed due to ferromagnetic 
spin fluctuations peaked at q = [BM]. Therefore it 



is natural to expect the origin of the triplet pairing in 
Sr2RuC>4 is also ferromagnetic spin fluctuation |,§. This 
assumption has been believed to be justified by NMR 
experiments j7|-[|. However the recent neutron scatter- 
ing experiment has shown that there exists a significant 
peak near q = (±27r/3, ±27t/3) and no sizable ferromag- 
netic spin fluctuation |10|. Thus it is difficult to assume 
that the spin fluctuation near q plays no role in the 
Cooper pairing in Sr2Ru04. (In the following discus- 
sion we call this fluctuation as antiferromagnetic (AF) 
spin fluctuation, for simplicity.) However this AF fluctu- 
ation leads to the singlet superconductivity rather than 
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the triplet superconductivity as expected in analogy to 
high-T c cuprates g. 

In this paper we propose a mechanism which gives the 
triplet pairing even if the spin fluctuation is AF. We find 
that the characteristic features of Sr2Ru04 are twofold: 
One is the anisotropy of the spin fluctuation found in 
NMR experiments [|]|]] , and the other is a nesting prop- 
erty with momentum q of the two-dimensional Fermi 
surface. We show that these two features explain the 
pairing in Sr2Ru04. 

In addition to the competition between singlet and 
triplet pairing, the direction of the d vector, which is the 
order parameter of triplet superconductivity, is another 
interesting problem. We show that the anisotropy of the 
spin fluctuation also explains the experimental fact that 
the d vector is parallel to the z-direction 0] . First we ex- 
tend the RPA formulation to the case of anisotropic spin 
fluctuation. Using the obtained effective interactions, we 
investigate the most stable pairing based on the weak- 
coupling gap equations. When the spin fluctuation is 
isotropic, the so-called d 2 .2_ y 2-wave pairing is the most 
stable. However when the anisotropy is increased, the 
state corresponding to z(sink x ± isiak y ), which is the 
prime candidate of Sr2Ru04, becomes the most stable. 

For the 7 band which is one of the three bands in 
Sr2Ru04 we assume a two-dimensional effective 

Hamiltonian 
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where Ck a is the annihilation operator of an electron with 
momentum fc and spin a. We consider only the on-site 
Coulomb repulsion, /, as in the previous studies of spin- 
fluctuation mechanism. Among the three bands, we con- 
sider the 7 band consisting of the antibonding band of 
Ru 4d a .j / and O orbitals in this paper, because it 
has the largest density of states at Fermi energy and the 
superconductivity is considered to be realized predomi- 
nantly in the 7 band 0| . Although the spin fluctuation 
near q is understood from the nesting effect of a and (3 
bands ||, we assume that the wave-number dependence 
of spin fluctuation is common in the three bands due to 
some interactions, such as spin-orbit couplings and/or 
Hund couplings. In the following calculation, we use the 
7 band. However the same discussions can be also ap- 
plied to the a and (3 bands. 

The anisotropy of spin fluctuation observed experi- 
mentally is implicitly included in the two-body Hamil- 
tonian, Hq. Our purpose is not to investigate the origin 
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of anisotropy in details but to examine the role of the 
anisotropy to Cooper pairing. Therefore we introduce a 
phenomenological parameter a by 



X(+- o)(<?) = "X(TT, o)(Q), 



(2) 



where X(TT- o)(<z) (X(h , o)(o)) is the unperturbed static 

susceptibility of z axis (xy plane), which originates from 
Hq. The parameter a represents the anisotropy of spin 
fluctuation and we take a < 1 since NMR experiments 
show that x(xx) < X(zz) §,§• 

Using this one-band model, we discuss the effective in- 
teractions between Cooper pairs due to spin fluctuations. 
Summation of bubble and ladder diagrams (i.e., RPA ap- 
proximation) gives 
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with 



V h . (k - k') = 

Vb.e(fc - k') = 



I (I/N) Xm , 0) (k-k') 

Ni-(i/Nr x 2 w , 0) (k-k'y 

I (I/NYxl ]n) {k-k') 
Nl-(I/Nf X 2 m>0) (k-k'y 



V l3A (k-k) = — — — -± '—- — -j-. 4 

N 1 - [I/N)x(+- I o){k-k) 

Here Vb.o (14>.e) comes from the summation of diagrams 
with odd (even) number of bubbles, and I4ad from the 
ladder diagrams. It is apparent that Vb.o is between the 
electrons with equal spins while Vt>. c and Mad are between 
those with the opposite spins. 

It is straightforward to derive the gap equations in the 
anisotropic case, using the method developed by Leggett 
Jl3[ . First we introduce the operators 



4^ = y^i^ss'C-ksCks' 



(5) 



4 = ^2(v2Va)ss>C- k sCks', for a = 1, 2, 3, (6) 



where o~ a (a = 1,2,3) are Pauli matrices. The opera- 
tor {t^) corresponds to spin singlet (triplet) Cooper 
pairs. In terms of these operators, the effective interac- 
tion (||) can be rewritten as 



kk' 

+iEE^ ) (fe-fe')4° )t 4' ) : (7) 



where 

V sin (k - fe') = 2[V b . c (fe - k') + Mad(fc - fe')]> 

V&\k fe') = V t [f(k - fe') = -2V h . (k - k'), 

V^{k - fe') ee 2[V h . c (k - fe') - V lad (fc - fe')]. (8) 

Since Sr2Ru04 has a long coherence length in ab plane, 
£ a b ~ 660 A Ej ], we use mean- field approximation to 
-f^int ■ We restrict the discussion to unitary states because 
it is unrealistic to assume non-unitary states in S^RuCU 
|l5| . Requiring that there is no coexistence of singlet and 
triplet pairs, we obtain the gap equations 

A(fc) = - £ V sin (k - fc')A(fc')e(£ sin (fc')), 

fe' 

dW(fe) = -^2v£\k - fc')d( Q )(fc')e( J B tri (fc')), (9) 

fe' 

where Q(E) ee ^tanh^, £ s 2 in (fc) =£ + A(fc)A*(fc), 
and £ t 2 ri (fc) = Cfc + d(fc) • d*(fc) with f fc = e fc - /x. 
The singlet and triplet order parameters are defined 

as A(fc) - -2Efe'Kin(fe - fe')(4' ) > and rf(a) ( fc ) = 



~5 Sfe' ^tri — ^')(4' )' respectively. Here d(fc) is the 
so-called d- vector for the triplet superconductivity. 
In the system with the rotational symmetry in spin 

space, X(TT- o)(q) = X(h , o)(*z) i s satisfied and thus the 

relation Vb.o + Vb. c = Viad holds. In this case, it is easy 



fefe' a=l 



to see V&>(k - k') = V^{k - fe') = V&\k - k')_ 

On the other hand, the gap equation in Eq. (H) for 
the triplet pairing becomes dependent on the direction 
of the d vector in the anisotropic case. It means that d 
vector has some preferred direction if the triplet pairs are 
formed by anisotropic spin fluctuations. This is naturally 
understood because the d vector is orthogonal to the spin 
direction of triplet Cooper pairs fl6| . For the present 
case with %(+-, o)(<z) < X(TT, o)0z) (i- e -> a < 1) which 
is applied to the Sr 2 Ru04, we can see from Eq. (f|) that 
Vi a d{k — fe') is suppressed and the effective interaction 
V^.j (fe — fe') approaches Vsi n (fc — fe'). Consequently the 
triplet superconductivity with d^(k) (i.e., d || z) can be 
stabilized even due to the AF spin fluctuations. 

In order to determine the symmetry of the supercon- 
ducting order parameter, we have to take account of their 
sign change along the Fermi surface. For the high-T c su- 
perconductors, the AF spin fluctuation with momentum 
(it, 7t) stabilizes the singlet d a ,2_ J( 2-wave superconductiv- 
ity. In that case, the singlet order parameters A(fe') with 
fe' = (tt, 0) and A(fc) with fe = (0,7r) have the opposite 
sign, so that the gap equation in (|9|) is satisfied with 
F sin (7r,7r) > 0. 

For Sr2Ru04 we consider that a kind of nesting prop- 
erty of the Fermi surface plays an important role. This 
is the second point of our mechanism. Figure 1 shows 
a schematic Fermi surface for the 7 band. Since the AF 
fluctuation in Sr 2 Ru04 has momentum q , the Fermi sur- 
face is also shifted by (2n/3, 2ir/3) in Fig. 1. It is appar- 
ent that some part of the shifted Fermi surface overlaps 
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with the original Fermi surface with modulo 2ti. In anal- 
ogy to the case of high-T c superconductivity, if the super- 
conducting order parameters have the opposite sign on 
these overlapping portions of the Fermi surface, the gap 
equation is satisfied with (2-7r/3, 27r/3) > 0. From 
Fig. 1, it is natural to consider the p-wave pairing in- 
stead of the singlet d x 2_j,2-wave pairing. 
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FIG. 1. A schematic Fermi surface for the 7 band of 
S^RuCU. In order to show the nesting property, the Fermi 
surface shifted by the AF wave number, go, is also shown 
by the solid line. The thin dashed lines indicates the Fermi 
surfaces in the extended Brillouin zone. 

In order to clarify this point quantitatively, we com- 
pare various kinds of anisotropic superconductivity using 
the effective interaction and the simplified Fermi surface. 
Near the transition temperature T c , we rewrite the gap 
equations as 



0(fc) = -J2 V*(k - k'Mk')-±- tanh ^ 



(10) 



where <f>(k) represents one of the order parameters A(fc) 
or S a \k), and is determined from Eqs. (@) depending 
on <j). In the weak coupling approximation, T c is obtained 
as 



ksT c — 1.13hvFk c exp 



1 



mm)) 



FS 



(11) 



where vf, k c and N(0) are the Fermi velocity, cut-off of 
the wave number, and the density of states at the Fermi 
energy, respectively. ((V^fs means the average over 
the Fermi surface, 



FS 



— JFS 



f FS dk J FS dk'V<p(k - k')(j)(k)(j)(k') 



[J FS dk']J FS dk^(k) 



(12) 



We identify that the order parameter which gives the 
largest N(0)((V^,))fs is realized. 

For Sr2Ru04 we choose order parameters 4>(k) as fol- 
lows 



61 (fc) = cos k x + cos k y , 

6 2 (fc) = cos k x — cos k y , 

63 (fc ) = sin k x sin k y , 

64 (fc) = sinfc x , (d _L z), 
65 (fc) = sin fc^ (d || z), 



(13) 



where <pi ~ (f> 3 correspond to singlet pairings, and 04, 05 
to triplet pairings, respectively. The most probable can- 
didate for Sr2Ru04 is z(sink x ± ismk y ) which is equiv- 
alent to </>5 just below T c , because the gap equation ( |To| ) 
for sin k x ± i sin k y is exactly same as that for 05 . If 
A r (O)((V0 5 ))i?s is the largest, we expect that the order 
parameter d(k) = z(sink x ±ismk y ) is realized, because 
near zero temperature it acquires a larger energy gap 
than 05. 

To emphasize the characteristic feature of the nesting, 
we assume the simplified Fermi surface as shown in Fig. 
1. For the q dependence of X(TT, o)(?) with a maximum 
at qo , we use the susceptibility obtained in the LDA cal- 
culation ||], and fix ,5(0) = 0.8 with S(q) = jjXtfh o)(<?)- 
We regard S(q ) as a phenomenological parameter. 

Figure 2 shows the a dependence of N(0)({V,f, n ))FS 
(n = 1 ~ 5) for S(q ) = 0.95. We examined various 
choices of S(q Q ) from 0.90 to 0.99 to find that the results 
do not change qualitatively. When the anisotropy is weak 
(a ~ l) ; the singlet d^ 2 ^2 -wave superconductivity, 02, 
is stabilized. On the other hand, when a is small, the 
order parameter 05 is stabilized which is consistent with 
experiments. 
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FIG. 2. The dependence of the anisotropy parameter, a, of 
iV(0)<(V^ n » FS (n = 1 ~ 5) for S(q ) = 0.95. 

The phase diagram as a function of a and S(q ) is 
determined by examining various values of S(q Q ). Be- 
cause it is unphysical to assume that S(q ) is very close 
to 1, we show the results up to S(q Q ) = 0.99 in Fig. 3. 
When the spin fluctuation is isotropic (i.e., a — 1), the 
singlet d a; 2_, y 2-wave superconductivity, cosk x — cosfcj,, is 
the most stable pairing. This is consistent with the previ- 
ous studies |q| . However we find a fairly large parameter 
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region where the state corresponding to z(sin k x ±i sin k y ) 
is realized. 




S(q ) 

FIG. 3. Phase diagram as a function of the anisotropy 
parameter, a, and S(qo) which is the maximum of 
S(g) = (I/N) X (n, o)(g). with q = (±2tt/3, ±2tt/3). 

Finally we discuss the competition between the singlet 
cos k x — cos k y pairing and the triplet z(sin k x ±i sin k y ) 
pairing in terms of the effective interaction and the nest- 
ing property. From the explicit form of V^ n for n = 2 
and 5, we can see that a relation V& 2 > Va, 5 is satisfied. 
Therefore if we consider only the magnitude of the effec- 
tive interaction, the singlet pairing is favorable. However 
the nesting property favors the triplet pairing. Let us 
assume that V^ n is enhanced very strongly by the AF 
fluctuation and approximated as 

(q) = jjAj(q x ± 27r/3)% a ± 2tt/3), (14) 

with S being the S function with modulo 2n. Using this 
approximated form of V ( j >n (n = 2 and 5), we obtain 

n (0)((V^))fs = [-2.79 x 1(T 2 <S(0) + 4.91 x lO" 2 ]^, 
N(0){(V^}) FS = [4.24 x 1CT 2 (5(0) + 5.06 x lCr 2 ]A 5 . 

This estimation shows that the z(sin k x ±i sin k y ) pairing 
utilizes the peak of X(ft, o)(o) a ^ Qo more effectively than 
cos k x — cos k y pairing does. Therefore, even if A 2 > A 5 , 
the triplet pair can be stabilized. 

In determining the phase diagram in Fig. 3, we have 
assumed simple functional forms of the order parameters, 
<j) n (k). For the detailed calculations, it will be necessary 
to optimize the fc-dependence of <f> n (k). However the 
global feature of the phase diagram will not change. 

In summary, we have generalized the RPA formulation 
of the effective interaction due to the spin fluctuations 
and derived gap equations including the anisotropic case. 
We have shown that the state corresponding to z(sin k x ± 
i sin k y ) becomes the most stable even if the AF spin 
fluctuation is dominant, when the anisotropy is strong 



enough and the nesting property of the Fermi surface 
is present. Although the nesting for the actual Fermi 
surface will be weaker than what we assumed here, it 
is reasonable to think that our mechanism is the most 
promising one as far as the AF fluctuation is dominant. 

In this paper we have investigated the pairing in the 7 
band. However it is straightforward to consider the other 
bands (a and /3 bands) in S^RuCi. Since the nesting 
property will be comparable or even stronger for these 
bands than for the 7 band, we expect the same mecha- 
nism for triplet superconductivity for a and j3 bands even 
if the 7 band does not have the peak near q . 

It is reported that Sr 2 Ru04 has exotic property called 
as 3K phase Jl7| when Ru metal is embedded in the single 
crystal. We speculate that the enhancement of T c is due 
to the increase of the anisotropy (i.e., decrease of a) near 
the interface region between Sr2Ru04 and Ru metal. We 
consider that to investigate the origin of the anisotropy 
is very important both for understanding the supercon- 
ductivity and for finding the new exotic phenomena. 
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thanks to H. Namaizawa for useful instructions. 
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